Dynamics of Fluxon Lattice in Two Coupled Josephson Junctions 
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We study theoretically the dynamics of a fluxon lattice (FL) in two coupled Josephson junctions. 
We show that when the velocity of the moving FL exceeds certain values (V a ,b), sharp resonances 
arise in the system which are related to the excitation of the optical and acoustic collective modes. 
In the intearval(Vo) Vb) a reconstruction of the FL occurs. We also establish that one can excite 
localized nonlinear distortions (dislocations) which may propagate through the FL and carry an 
arbitrary magnetic flux. 



In recent years many papers have addressed the dy- 
namics of the fluxon lattice in layered superconductors 
and in particular, in high T c superconductors. In the 
absence of the magnetic field, the spectrum of collective 
oscillations in the long wave limit has a threshold fre- 
quency which coincides with the Josephson frequency wj 
- 0. If a magnetic field is applied, the vortex lattice 
is formed in a superconductor, and the spectrum of col- 
lective modes is changed. In particular, an acoustic-like 
mode related to the vortex lattice oscillations arises. The 
spectrum of the fluxon lattice (FL) arising in a parallel 
magnetic field was calculated in Refs. || - ||. 

The collective modes may be excited by an external ac 
field and also by a dc current (jdc) across the layers. In 
the latter case, the collective modes are excited if the ve- 
locity of the moving FL coincides with the phase velocity 
of the collective modes. This effect has been studied in 
detail for the case of a long Josephson tunnel junction ||] 
and was studied in a theoretical paper |l(| recently for 
the case of layered superconductors. 

In the present paper we consider two coupled, long 
Josephson junctions and study the excitations of small 
amplitude collective modes as well as the excitations of 
nonlinear perturbations (dislocations) of the FL in such a 
system by a dc current. Comparative simplicity of equa- 
tions governing the FL dynamics in this system, allows 
one to analyse effects arising in this system in detail and 
to understand the behaviour of more complicated struc- 
tures, like layered superconductors. It will be shown in 
particular that nonlinear excitations (dislocations) may 
arise in the system by the dc current and that these ex- 
ciations can carry an arbitrary magnetic flux (larger or 
smaller than the magnetic flux quantum <E>o)- We note 
that the system under consideration was analysed in the 
absence of a magnetic field in Refs. [[llj - 

I Model and Basic Equations 

Let us consider the system shown in fig.[l| We as- 
sume that different currents may be passed through the 
junctions 1 and 2, i.e. a current through the middle su- 
perconductor can be driven independently from currents 
through the outer superconducting electrodes. For sim- 
plicity we assume that the junctions are identical, i.e. 
they have equal critical currents etc. Equations describ- 
ing the dynamics of the coupled Josephson junctions were 



obtained in a number of works (§ - g, Q - Q. The 
magnitudes of the magnetic field in junctions 1 and 2 are 
related to the phase difference (p% t 2 through the expres- 
sion 



#1,2 = (l/2)d x [<^i,2 + 7^2, l] 



(1) 



Here i?i,2 are the dimensionless magnitudes of the 
magnetic field in the junctions 1 and 2. We choose the 
quantity Hq = $>o/2ir\ 2 as a unit of measurement of the 
magnetic field and the London penetration depth as a 
length unit. $o is the magnetic flux quantum, tpi^ is the 
phase difference in junctions 1 and 2 and 7 = exp[-2d], 
where 2d is the thickness of the middle electrode (in units 
of A). We assume that the characteristic scale of spatial 
variations of ip± t 2 is much greater than unity, (i.e. than A) 
and also that the thickness of the outer superconducting 
layers are greater than A. 

The current through the junctions can be expressed 
through the corresponding component of (V x H) and 
be related to the quasiparticle and Josephson currents. 
We obtain 



2ljd x H 



1,2 



= (df t + ad t ) fi,2 + sin(ip lt2 ) - 771,2 (2) 

Here lj = {cHq/SttJcX) 1 ^ 2 is the dimensionless Joseph- 
son penetration length, j c is the critical Josephson cur- 
rent density, a = h/2ep qp j c to is the damping constant, 
p qp is the junction resistivity due to quasiparticle tunnel- 
ing. Time is measured in units to = y/^C/2ej c , where C 
is the junction capacitance (per unit area). The constants 
771,2 are dimensionless currents (in units of j c ) through 
the junctions 1 and 2. It can easily be shown that the 
magnetic flux in the system equals 



dxd x (ipi + (f 2) 



(3) 



Substituting for H12 in (Eh from 
of two coupled equations for ipi 2 



1), we obtain a set 



Ijdxx Vpi,2 + 7^2,1] = (dtt + ad t ) ^1,2 + szn(^i, 2 ) - 771,2 

(4) 

We introduce the new functions tp± — (l/2)(yi ± ip%). 
Summing and subtracting eqs.(^), we obtain the new 
equations for tp±. 
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l±d xx ip± = {d tt + ad t )(fi± + sin((p±)cos(ip T ) - rj± (5) 

Here l\ = ± 7) and r]± = (771 ± m )/2. Eqs.(§) de- 
scribe the dynamics of two coupled Josephson junctions. 
We use them for studying the FL. Let us assume that 
a magnetic field parallel to the planes of the Josephson 
junctions is applied and a dense FL arises in the junc- 
tions 1 and 2. In the stationary state (and sufficiently 
high magnetic fields) the solution for ipi^ can be easily 
found from eqs.(|): tpQ = Hx ± (tt/2) + 4y\ Where 

w ^{l-TL)~ 2 cos{TLx). This solution is valid pro- 
vided 



Here AA± = \D_\~ 2 ± \D+\~ 2 , 4B = bJD. 



C 1 <7i<2 7 r 



(6) 



Then the right hand side of this condition means that 
eq.(Q) is valid, i.e. the characteristic scale of ipi^ varia- 
tion along the x-axis is greater than the London penetra- 
tion depth. The left hand side of this condition means 
that the FL is dense and spatial oscillations of (pi t 2 are 
small. The field in the junctions TL is related to the ex- 
ternal field by TL — 2H e / (1 -I- 7). The expression for 
given above describes two fluxon chains each of which is 
shifted by a half period with respect to each other. 

II Dynamics of the Fluxon Lattice 

Let us now consider solutions describing the motion 
of the FL driven by the dc currents 771,2 (the currents 
7/1 1 2 may differ from each other). We seek the solution of 
eqs.(|^) in the form of a traveling wave assuming that the 
junctions are long enough and neglecting reflected waves 



Wl,2 = Ux -Vt + 7^1.2 + 01,2 + 



9(0) 
'1.2 



(7) 



Where ipip is the rapidly oscillating part of y>x,2 m 
space with a period of 2n/TL. #1,2 is the slowly varying 
part and 9{\ = (for 1) and tt (for 2). A similar repre- 
sentation was used in Ej] - Jl6| for finding the shape of 
the supersolitons (dislocations) in the FL and in M for 
finding the spectrum of acoustic-like oscillations of the 
FL. For the rapidly oscillating part ip± = (-01 ± V'2)/2, 
we have from eqs.(|5|) 

= - tS^F i b + C0S ( Y + e +) + aVsin(Y + 0+)} 



\D+i 
cos(0_) 



{-bsin(Y + 0+) + aVcos{Y + 0+)} 



(8) 



where Y = Hx - Vt and 0± = (0i ± 2 )/2, b± = V£ 
- V 2 , D± = b± + iaV and V± = l±H. Assuming that 
ij)± are small (i.e. j <C 1), and expanding sn\(ip±) in 
powers of ip±, we obtain from eq.(^) the equations for 
the slowly varying part 0± in the main approximation 

l 2 + d 2 xx e+ = (dl + ad t )9+ + aV A-cos(26 L) 
+aVA + +aV -f]+ 
l 2 _d 2 x J- = (d 2 t + ad t )9- + Bsin(29-) - rj- (9) 



b+\D + 



Eqs.(^|) describe the dynamics of the FL. The 
phase 0+ is a local displacement of the FL as a whole and 
the phase 0_ determines a relative displacement of two 
fluxon chains. 

II a. Collective modes 

Consider the stationary case when the currents 77^ are 
absent and the FL is motionless (V = 0). Lineariz- 
ing eqs.(H), we obtain for the spectrum of the collective 
modes (for simplicity we neglect the damping) 



= K I 



+ > 



= U)n+ K 2 f 



(10) 



Here cu 2 — r y/[(l — r ) 2 )l 2 'H 2 ] is the threshold frequency 
for the optical branch. It decreases with increasing TL. 
The first expression in eq.(pTj) describes the acoustic 
branch of the FL oscillations. Similar modes exist in 
layered superconductors || - ||. These modes are inde- 
pendent from each other. In the acoustic (optical) mode 
the phase 0+ (0_) are perturbed. Consider now excita- 
tions of these modes by the moving FL. 

II b. Moving fluxon lattice 

Let us assume that a current 771 = 772 flows through 
both junctions (then 77+ = 771 = 77 and = 0). Then the 
stable solution of eq.(|) is: 20_ = for V < V a , V > V b 
and 20_ = tt for V a < V < Vb- Here V a ^ are the roots 
of the equation B(V) = 0. The solution 0_ = 0, (tt) 
corresponds to a positive, (negative) value of B. For the 
case of a small damping a (i.e. a <C V±) we have V 2 b = 
Vl[lTa 2 /{V 2 -V 2 )] 

If the velocity of the moving FL exceeds (V a /H), a 
reconstruction of the moving FL occurs (see fig.||). If the 
velocity of the FL increases further and exceeds (Vb/H), 
the initial triangular form of the FL is restored. 

The form of the current-voltage characteristics ?7(V) 
may be easily found from eq.(^) averaged in space. We 
obtain 



//-M- = 2nr<; j^-j-2 



v < v a , v > v b 
v a < v < v b 



(ii) 



There are two peaks in the ?7(V) characteristics, one 
of them corresponds to the velocity of the optical mode 
(at large k) and the second corresponds to the velocity 
of the acoustic mode. In addition there are two jumps 
at voltage V = V a and V = Vb corresponding to the FL 
reconstruction. If the damping constant a is small, then 
V a ,b ~ V±. In fig.||we present the form of the r](V) curve 
assuming for simplicity that a does not depend on V. 

II c. Creation of dislocations 

Let us assume that a bias current flows in the mid- 
dle superconducting electrode (?7_ ^ 0). One can see 
that the second equation of (^]) is similar to an equation 
describing a single Josephson junction. This equation 
describes nonlinear distortions in the FL. These distor- 
tions (dislocations) were analysed in Jl4[ where it was 
shown that they arise as kinks in the FL in two, slightly 
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different coupled Josephson junctions. Similar distor- 
tions (supersolitons) may arise in a single long Joseph- 
son junction whose parameters are modulated in space 
p5| . The characteristic length of dislocations l_/B de- 
pends on both the magnetic field and the applied voltage 
(or current). Dislocations are created at r\- > B(V,Tt) 
and can propagate through the system in an interval of 
?7_ above and below B (Fiske steps If the damp- 

ing is small, the dislocation has the well known fluxon 
form 2#_ = 4tg~ 1 exp(£), where £ = (x - ut) /Id and Id 
= l-/\f~B. The velocity u is related to via the well 



known formula |16| which in the slow velocity limit is 
reduced to u = (ir/4:)r]^ld/a. According to eq.((|), a dis- 
location arising in a moving FL causes a perturbation of 
the phase 9 + . Substituting the expression for into (||), 
we obtain for 8 + 



cosh £ 



(12) 



Where (3 = aul d /[l\~u 2 ] ands = 2aFA_ ■l 2 d /[l\-u 2 ]. 
If the FL velocity V lies in the intervals (0,V a ), (14, oo), 
we obtain from eq.(pl) rj + = aV{A + + A-). 



From eq.(12) we obtain 



d^icosh £1 • e 



/3(«i-«) 



(13) 



The spatial dependence of the dislocation 6* + (£) is 
shown in fig. || for different values of (3. The expression for 
0+{£) (13) is valid provided that the characteristic size of 
the dislocation (3~ l is less than the junction length L. Let 
us calculate the magnetic flux carried by a dislocation in 
a moving FL. Substituting expression (p^) into eq.(^), we 
have for magnetic flux in the system 



$ = 2(HL) + 2 J dx d x 0-\ 



(14) 



The first term in eq.(|l4|) is the magnetic flux in the 
system in the absence of a dislocation. The second term 
is the magnetic flux $d carried by a dislocation. With 
the help of eq.((l^) we obtain for ^ 



8VA- 
u{V) 



(15) 



One can see that the flux &d depends both on the ve- 
locity of the dislocation and the velocity of the FL as a 
whole, turning to zero at V = 0. Therefore, a dislocation 
in the FL is a localized distortion which can move un- 
der the action of an external force (the difference of the 
currents ?7_) and carry an arbitrary magnetic flux, the 
magnitude of which is determined by currents r) + and 

III Conclusions 

We analysed the dynamics of a dense FL in a system 
of two coupled Josephson junctions. Acoustic and opti- 
cal collective modes may propagate in the FL. If the FL 



is moving in the presence of a dc current, a resonance 
excitation of the modes takes place when the FL veloc- 
ity coincides with the limiting velocity V- of the optical 
mode, or with the velocity of the acoustic mode V+ . In 
the interval (V-,V + ) the FL is reconstructed. It was also 
shown that if the currents r\\ and 772 through the junctions 
are different, localized distortions (dislocations) may be 
created in the FL. They carry an arbitrary magnetic flux 
and may lead to a non- Josephson generation. 
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FIG. 1. The structure under consideration (schematically). 
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FIG. 2. The deviation of the current from Ohm's Law 
(8-q — rj — aV) due to excitation of the collective modes vs 
voltage. (We used a = 0.5) The positions of fluxons (crosses) 
in both junctions are shown inset for V > V a ,Vb < V (a) and 
V a <V <V b (b). 
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FIG. 3. The spatial dependence of the magnetic field (9' + 
— d^6+) in the dislocation for different values of f3: (5 — 0.2 
(1), 0.2 (2), 1.0 (3). 
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